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ABSTRACT 

Outflows and jets are believed to play a crucial role in determining the mass of the central protostar 
and its planet forming disk by virtue of their ability to transport energy, mass, and momentum of 
the surrounding material, and thus terminate the infall stage in star and disk formation. In some 
protostellar objects both wide angle outflows and collimated jets are seen, while in others only one is 
observed. Spitzer provides unprecedented sensitivity in the infrared to study both the jet and outflow 
features. Here, we use HiRes deconvolution to improve the visualization of spatial morphology by 
enhancing resolution (to sub-arcsecond levels in the IRAC bands) and removing the contaminating 
sidelobes from bright sources. We apply this approach to study the jet and outflow features in Cep 
E a young, energetic Class protostar. In the reprocessed images we detect: (i) wide angle outflow 
seen in scattered light; (ii) morphological details on at least 29 jet driven bow shocks and jet heads or 
knots; (iii) three compact features in 24 fim continuum image as atomic/ionic line emission coincident 
with the jet heads; and, (iv) a flattened ~35'' size protostellar envelope seen against the interstellar 
background PAH emission as an absorption band across the protostar at 8 fim. By separating the 
protostellar photospheric scattered emission in the wide angle cavity from the jet emission we show 
that we can study directly the scattered light spectrum. We present the H2 emission line spectra, as 
observed in all IRAC bands, for 29 knots in the jets and bowshocks and use them in the IRAC color 
- color space as a diagnostic of the thermal gas in the shocks driven by the jets. The data presented 
here will enable detailed modeling of the individual shocks retracing the history of the episodic jet 
activity and the associated accretion on to the protostar. The Spitzer data analysis presented here 
shows the richness of its archive as a resource to study the jet/outflow features in H2 and scattered 
light in a large homogeneous sample. 

Subject headings: infrared: general — star formation: protostar — ISM: jets and outflows 



1. INTRODUCTION 

Star formation begins with the collapse of a dense in- 
terstellar cloud core to a protostar surrounded by a disk 
of gas and dust. The envelope of the cloud core falls 
inwards to feed further growth of the protostar and its 
accreting disk. At some point during the accretion phase, 
protostellar jets and winds originate close to the surface 
of the forming star and interact with the dense enve- 
lope surrounding the protostar-disk system (e.g. Konigl 
& Pudritz 2000; Shu et al. 2000). These outflows and 
jets are believed to play a crucial role in determining the 
mass of the central protostar and its planet forming disk 
by virtue of their ability to transport energy, mass, and 
momentum of the surrounding material, and thus ter- 
minate the infall stage in star and disk formation. It is 
now well recognized that jets and outflows are essential 
and inherent to the star formation process. These flows 
are broadly classifled into two types: molecular outflows 
traced mainly with CO emission, and jets observed by 
optical emission (cf. recent reviews by Arce et al. 2007, 
Pudritz et al. 2007; Bally 2007). Observations indicate 
that around some protostars, high- velocity jets with a 
narrow opening angle are enclosed by a low- velocity out- 
flow with a wide opening angle (e.g. Mundt & Fried 1983; 
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Velusamy et al. 2007) however, in others only one com- 
ponent is observed. Poorly collimated flows can be due to 
extreme precession of the jet (Shepherd et al. 2000) and 
indistinguishable from wide angle outflows. The wide 
angle outflows, which are also observed in the scattered 
light, are often missed in the optical or near-IR due to 
instument insensitivity to faint emission and confusion 
from the bright protostellar emissions. Spitzer provides 
unprecedented sensitivity in the infrared to detect the 
jet/outflow features. The IRAC, IRS, and MIPS obser- 
vations can provide new insight into the structure, mor- 
phology, and physical and chemical characteristics of the 
outflow sources as demonstrated in the observations of 
HH46/47 by Noriega-Crespo et al. (2004a) & Velusamy 
et al. (2007), Cep E by Noriega-Crespo et al. (2004b), 
L1512F by Bourke et al. (2006), L1448 by Tobin et al. 
(2007) & Dionatos et al. (2009), L1251 by Lee et al. 
(2010), and HH211 by Dionatos et al. (2010). Using the 
IRAC colors alone, Ybarra & Lada (2009) demonstrated 
the feasibility to study the pure rotational H2 line emis- 
sions in the shocks without the need for spectroscopic 
data. Velusamy et al. (2007) demonstrated that resolu- 
tion enhanced re-processing of Spitzer images in the case 
of HH46/47 brought out more clearly the morphologies 
of: (a) wide-angle outflow cavities in the IRAC 3.6 & 4.5 
/im images by scattered starlight, (b) jets and bow shocks 
in H2 molecular emission within the IRAC bands, and 
(c) the hottest atomic/ionic gas in the jet head charac- 
terized by [Fe II] and [S I] line emissions within the MIPS 
24 /im band. 
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In this paper we present new results on Cep E obtained 
by reprocessing the Spitzer IRAC and MIPS images with 
the HiRes deconvolution (Backus et al. 2005). Cep E is a 
remarkably bright outflow object, at a distance of 730 pc, 
with strong radiative power (~ 80 Lq), and appears to be 
an intermediate mass object which is expected to reach 
a final stellar mass of 3Mq (Froebrich et al. 2003). Cep 
E has been recently classified as a Class protostar by 
Noriega-Crespo et al. (2004b) based on the Spitzer EdLily 
Release Observations. It is among the youngest proto- 
stars with an age in the range of 20 - 400 thousand years, 
estimated for different evolutionary models (Froebrich, 
2005; and the references therein). In the near- and mid- 
IR its morphology is similar to that expected for a jet 
driven outfiow, where the leading bow shocks entrain and 
accelerate the surrounding molecular gas. The properties 
of this outflow have been well studied by Moro-Martm et 
al. (2001) and Smith et al. (2003). The southern bright 
bowshock, identified as HH 377 (Devine, Reiputh & Bally 
1997), is well studied in the optical wavelengths using the 
[O I] A 6300, [NI ] A 5200 and [S II] AA6717,6731 lines 
tracing the dissocitive shocks characteristics of HH flows 
(e.g. Ayala et al. 2000) and later in the near-IR in the 
[Fe II] 1.26 and 1.64 jam emission lines, utilizing the fact 
that [S II] and [Fe II] have very similar excitation energies 
(Nisini et al. 2002). Furthermore, Hartigan, Raymond & 
Hartmann 1987) interpreted the relative strength of [O I] 
and [N I] with respect to Hq, and the lack of [O II] and [O 
III] in terms of relatively low velocity (20-25 km s~"^) dis- 
sociative shocks and concluded that the south bowshock 
of Cep E is a low excitation HH object. In addition to 
dissociative shocks, the HH outflows are known to have 
also magnetic precursors (C-shocks) needed to explain 
excitation of the molecular H2 emission. Cep E lobes 
have been well imaged in H2 emission in IR: in v=l-0 
S(l) 2.122/im; v=2-l S(l) 2.248 /im; v=3-2 S(3) 2.202 
fim (Eisloffel et al. 1996; Ladd & Hodapp, 1997); v=0-0 
S(5) 6.91 jam and S(3) 9.66 /im (Noriega-Crespo 1998); 
v=l-0 S(l)1.21 /im. In the mid-IR a series of 6 pure 
rotational emission lines H2 0-0 S(2) to S(7) in the wave- 
length range of 12.28 /im to 5.51 jam in the Cep E lobes 
were first imaged by ISOCAM CVF (Moro-Martm et al. 
2001) and their excitation diagrams led to the use of H2 
line ratios as C-shock diagnostics. The emissions in the 
Spitzer images in all the IRAC bands are dominated by a 
wide range of H2 rotational transitions from 0-0 S(4) to 
0-0 S(13) and a few vibrational 1-0 transitions in the 3.6 
/im band. Thus Spitzer observations of Cep E enhance 
the rich data set already available in optical, near-IR and 
mm wavelengths. Furthermore, all the interpretations so 
far have focused on treating this object as having jets 
and bow shocks, but no wide angle outflow. Because of 
the limited spatial resolution even the CO outflow mor- 
phology has been interpreted as entrained low velocity 
molecular gas associated with the bow shock driven by 
the jets (Moro-Martm et al. 2001). Our enhanced image 
analysis show the simultaneous presence of wide angle 
outflow along with the bow shocks powered by the jets. 
In this paper we highlight the morphological features in 
the low velocity wide angle outflow and the high velocity 
jets and discuss a wide range of structures from compact 
knots to extended envelopes as seen in the reprocessed 
Spitzer IRAC and MIPS images. 



2. ANALYSIS 

Noriega-Crespo et al. (2004b) presented the mosaic 
images of Cep E in the Spitzer IRAC bands at 3.5 /im, 4.5 
/im and 8 /im in a three color representation and the mo- 
saic image of MIPS band at 24 /im. However, some of the 
structures were missed in these mosaic maps due to ei- 
ther lower spatial resolution or diffraction lobes. Spitzer 
provides unprecedented sensitivity in the infrared, but 
the spatial resolution is limited by the relatively small 
aperture (0.85 m) of the primary mirror. To maximize 
the scientific return of Spitzer images we use the HiRes 
deconvolution processing technique that makes the op- 
timal use of the spatial information in the observations. 
The algorithm, "HiRes" and its implementation has been 
discussed by Backus et al. (2005) and its performance on 
a variety of astrophysical sources observed by Spitzer is 
presented by Velusamy et al. (2008). The HiRes de- 
convolution algorithm is based on the Richardson-Lucy 
algorithm (Richardson 1972; Lucy 1974), and the Max- 
imum Correlation Method employed by Aumann et al. 
(1990) for IRAS data. As demonstrated by Velusamy 
et al. (2007; 2008) the HiRes deconvolution on Spitzer 
images retains a high fidelity preserving all the main fea- 
tures. HiRes deconvolution improves the visualization 
of spatial morphology by enhancing resolution and re- 
moving the contaminating sidelobes from bright sources. 
The benefits of HiRes include: (a) enhanced resolution 
- 0.6'' - 0.8'' for IRAC bands; - 1.8" and -7" for MIPS 
24 /im and 70 /im, respectively; (b) the ability to detect 
sources below the diffraction-limited confusion level; (c) 
the ability to separate blended sources, and thereby pro- 
vide guidance to point-source extraction procedures; (d) 
an improved ability to show the spatial morphology of 
resolved sources. We reprocessed the images at all the 
IRAC bands and the MIPS 24 and 70/im bands using all 
available map data in these bands in the Spitzer archives 
containing the prostellar object Cep E centered at {a 
= 23h03ml2.5s, S = -h61°42'30"[J2000]). We used the 
pipeline processed basic calibrated data (BCD) down- 
loaded from the Spitzer Science Center (SSC) archives. 
We constructed mosaic images at each band using the 
MOPEX software (provided by SSC). In addition, we 
also produced reprocessed images at each band by apply- 
ing the HiRes deconvolution on BCDs following the steps 
outlined by Velusamy et al. (2008). The mosaic and the 
HiRes images cover an angular size of 5' centered on the 
Cep E protostar. All the HiRes images shown in this 
paper were obtained after 50 iterations unless mentioned 
otherwise. 

The mosaic and the HiRes deconvolved images in the 
IRAC and MIPS bands, shown in Figures 1 and 2 re- 
spectively, bring out the merits of HiRes processing. In 
these figures the top panels show the reprocessed HiRes 
images and the corresponding observed (mosaic) images 
are shown in the lower panels. The mosaic and HiRes 
images have identical angular extent and scale. The im- 
ages are presented using a square-root intensity scale and 
color stretches. The color stretches used are similar in all 
images and an example of the colors representing the rel- 
ative intensities is shown by the wedge at the upper left. 
In the MIPS 24 /im mosaic image only the protostar is 
obvious. However, in the HiRes image at least three ad- 
ditional features along the atomic jet become identifiable 
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as sources distinct from the "side lobes" present in the 
mosaic image. These were not identified in previously 
published mosaic images. 

The MIPS images shown in Fig. 2 have small pointing 
differences (~ 1^') with respect to the IRAC images in 
Fig. 1. A comparison of the point sources in the MIPS 
24 /im image with their counterparts in the IRAC bands 
indicated that the pointing at 24 fim was off with respect 
to the IRAC bands by -0.07^ in RA and +1.1'' in Dec- 
lination. We shifted the MIPS 70 /im map position by 
matching the protostar position with the MIPS 24 /im 
image. In all subsequent displays and analysis we use 
the MIPS images corrected for their position shifts with 
respect to the IRAC images. 

3. RESULTS AND DISCUSSION 

The Spitzer IRAC and MIPS bands offer a unique re- 
source to study a wide range of components simultane- 
ously: protostars, protostellar disks, outflows, protostel- 
lar envelopes and cores. At the short wavelength IRAC 
bands the outflow cones are observable in scattered light 
from the protostar through the cavity created by the jets 
and outflows (cf. Tobin et al. 2007; Velusamy et al. 
2007). A signiflcant fraction of the emission in the IRAC 
bands is also considered to contain emission from the 
H2 rotational lines (e.g. Noriega- Crespo et al. 2004a,b; 
Neufeld & Yuan, 2008; Smith & Rosen, 2005) and there- 
fore is an excellent tracer of H2 emission in the protostel- 
lar jets. Recently Velusamy et al. (2007) have shown that 
molecular jets and molecular gas in the bow shocks are 
readily identifiable in the IRAC bands, while the hottest 
atomic/ionic gases in the bow shocks, are also identifiable 
in the MIPS 24 /xm band which covers a few atomic/ionic 
emission lines. The dust emission from the protostar in 
the MIPS bands is a good diagnostic of the circumstel- 
lar disks. Finally, an extended flattened infall envelope 
perpendicular to the outflow has been seen in the IRAC 
8 /im in absorption against the background interstellar 
PAH emission (Looney et al. 2007; Tobin et al. 2010). 
Spitzer' s coverage of a broad range of IR emissions in the 
IRAC and MIPS bands with high sensitivity and photo- 
metric stability along with the IRS data in some cases, 
provide suflicient data for a comprehensive modeling of 
the SED to derive the physical characteristics and the 
evolutionary stages of protostars (cf. Robitaille et al. 
2007; Forbrich et al. 2010). However, by applying the 
HiRes re-processing we can extract even more informa- 
tion from the Spitzer data. The sensitivity, the reso- 
lution enhancement and removal of the diffraction lobe 
confusion have led to visualizing and characterizing more 
clearly the following features in the Cep E protostar: 

1. very high dynamic range maps in which the proto- 
star is well resolved on a sub-arcsec scale from the 
surrounding jets and outflows 

2. wide-angle outflow cavities in the IRAC 3.6 and 4.5 
/im images identifled by the scattered photospheric 
starlight; Though they are visible in the Mosaic 
(standard PostBCD products), they are brought 
out more clearly in the HiRes images. 

3. jets and bow shocks in H2 molecular emission 
within the IRAC bands. 



4. hottest atomic/ionic gas in the jet head traced by 
the [Fe II] /[S I] line emissions within MIPS 24 /im 
continuum band. 

5. flattened infall envelope around the protostar seen 
in extinction at 8 /im. 

3.1. Protostar 

The protostar in Cep E is detected in all IRAC and 
MIPS bands. The IRAC HiRes images shown in Fig. 1 
have high dynamic range, but the brightest features are 
shown saturated as they are optimized to highlight the 
low surface brightness emission. In all bands the highest 
brightness is at the protostar location while the jets and 
the outflow features are at very low brightness, at levels 
of ~ 1% or less of the peak brightness. The color stretch 
used in Fig. 1 is designed to bring out these low sur- 
face brightness features and therefore, underrepresents 
the peak brightness at the protostar. In Figure 3 we show 
the HiRes intensity contour maps of the protostar in all 
bands. To bring out the huge intensity contrast between 
the brightest emission at the protostar and the weaker 
emission of the outflow-jet features plotted in Fig. 3, we 
show selected contours: the lowest level at 0.01% and the 
highest at 50% of the peak. The observed flux densities 
and sizes for the protostar are summarized in Table 1. In 
the IRAC wavelengths the emission from the protostar 
appears to be unresolved with sizes comparable to the 
HiRes beam. The protostar flux densities measured in 
the HiRes maps listed in Table 1 are consistent with the 
photometric estimates discussed below in section 3.4. 

Cep E is believed to be a double source separated by 
1.4'' (-1000AU) as shown by the 222 GHz millimeter in- 
terferometric observations of Moro-Martm et al. (2001). 
However, as seen in Fig. 3, the overlay of mm source 
positions on the IRAC 3.6 /im map does not show any 
evidence for a double source in the IR bands. The sub- 
arcsec resolution (~ 0.6'' - 0.8") in the HiRes images 
would have clearly resolved any double source if present. 
Though the protostar emission is barely resolved in the 
maps, there is a hint of slight elongation roughly along 
the separation between the mm-double sources. Our re- 
sults do not show any evidence for the presence of a sec- 
ond protostellar component. Instead, we suggest that 
the resolved mm double source structure may represent 
circumstellar cold dust within a few x 10^ AU. 

3.2. Wide Angle Outflows and Jets: IRAC Images 

It can be diflicult to observe the extended low surface 
brightness features associated with the outflow or jets 
as they may be severely confused by the high brightness 
of the protostar and the disk. Deconvolution at best re- 
moves, or at least minimizes, the diffraction effects of the 
brightest features, thus enabling improved visualization 
of the low surface features around them. Furthermore 
the resolution enhancement provides a sharper view of 
the outflow cavity walls, molecular jets, and bow shocks. 

3.2.1. Wide Angle Outflow Component 

The radiation from deeply embedded young protostars 
is completely opaque at short wavelengths and we can 
observe only the scattered light escaping along the bipo- 
lar cavities produced by their jets and outflows. The 
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scattering models (cf. Whitney et al. 2003a & b; To- 
bin et al. 2007) show that such outflow cavities can be 
traced in Spitzer images in IRAC 3.6 and 4.5 fim bands, 
as is seen more clearly in the HiRes deconvolved images 
(Velusamy et al. 2007; 2008). In Fig. 1 a wide-angle 
biconical outflow is clearly detected in the HiRes images 
at 3.6 and 4.5 jam. Though its presence is marginally ev- 
ident in the mosaic images, its morphological details are 
delineated more clearly only in the HiRes images. The 
wide angle outflow is clearly detected in the light from 
the protostar scattered by the outflow cavity. Scattered 
light is brightest at the 3.6 and 4.5 jam bands. As ex- 
pected for the photospheric light emitted by a protostar 
at T ~ 4000K (see section 3.4), the scattered light trac- 
ing the wide angle outflow component is relatively weak 
at 5.8 /im and nearly absent in the 8 fim image. In con- 
trast to this wide angle outflow traced by the scattered 
light from the protostar, the narrower and collimated jets 
and the associated bow shocks are detected by their H2 
emission in all the IRAC bands. All the four IRAC bands 
contain a few H2 spectral lines (see below for details). 
While bands 1 and 2 show both the wide angle outflow 
(in scattered light) and narrow jet features (in H2 emis- 
sion), the bands 3 and 4 are dominated by the latter. 
In the HiRes images the jet emission features show well 
resolved knots and bow shocks. In the north lobe alone 
there are at least 16 such features (see Fig. 5a). 

The HiRes processing brings out the morphology of the 
jet and outflow more clearly. In Fig. 4 we show a compar- 
ison between the 4.5 fim and 8 fim images by overlaying 
a schematic of the geometry of the wide angle outflow 
with its inner boundary matching the outer extent of the 
jet and bow shocks. The 8 /im emission delineates the 
narrow jet traced by the entrained H2 emission, while at 
4.5 /im, in addition to the jet, the wide angle outflow cav- 
ity is clearly visible. The scattered light seen in the 3.6 
/im (Fig. 1) and 4.5 jam (Figs. 1 & 4) delineates a sharp 
conical boundary for the outflow cavity. The projected 
opening angle of the cavity is ~ 85° in the 3.6 /im im- 
age which is signiflcantly large and places it among those 
with very wide outflows (Velusamy & Langer 1998; Arce 
& Sargent 2006) and those with simultaneous presence 
of high velocity jets (Velusamy et al. 2007; Torrelles et 
al. 2010 and references therein). The wide angle out- 
flow seen here in the scattered light closely resembles the 
^^CO(l-O) low velocity lobes (Moro-Martm et al. 2001) 
with the blue shifted SW lobe approaching the observer. 
However, the spatial resolution in the CO maps (8^^ x 
7^') cannot resolve the sharp edges of the wide angle out- 
flow cavity walls as seen more clearly in the scattered 
light images with sub-arcsec resolution presented here. 

Poorly collimated flows can be due to extreme pre- 
cession of the jet (e.g. Shepherd et al. 2000) and are 
indistinguishable from the apparent wide angle outflows. 
In the Spitzer image, shown in Fig. 4(a), there is a "sec- 
ond jet /outflow" nearly perpendicular to the main north- 
south outflow as indicated by the east- west ellipse. The 
presence of this secondary jet was flrst observed, earlier 
in the H2 2.12 /im image (Eisloffel et al. 1996). The 
two brightest knots in the 4.5 /im image were clearly de- 
tected in their H2 2.12/im image while the eastern most 
knot is outside the 2.12 /im image. In a later H2 2.12 /im 
image (Ladd & Hodapp, 1997) all three 4.5 /im knots 
are detected. The presence of two outflows nearly per- 



pendicular to each other suggests the driving source in 
Cep E is a binary with a separation 4-20 AU (0.005^^ - 
0.03'0, and disk radius 110 AU (Terquem et al. 1999). A 
binary systems enhances the effects of precession of both 
the primary and secondary jets. Eisloffel et al. interpret 
the wiggles in the loci of the H2 2.12 jam knots as due 
to precession with a period of 400 yr and a pitch angle 
4°. Binaries with multiple outflows may be common; for 
example: HH111/HH121 (cf. Gredel & Reipurth, 1993) 
andL1551 (cf. Rodriguez et al. 1998). Therefore we can- 
not rule out that the observed wide angle cavity in Cep 
E is the result of precessing jets. Indeed the combined 
effect of the precessing jet along with the geometry of the 
"infall envelope" (discussed in section 3.5) is an alternate 
scenario as the cause for the observed wide angle cavity. 
In the absence of direct observational evidence on the 
outflow motions using high spatial resolution data of CO 
or other molecular tracers we cannot identify with cer- 
tainty the wide angle cavity detected in scattered light 
as wide angle outflow. 

In the jet driven outflow models the outflow emission 
near the protostar may be interpreted as the wings of the 
bowshocks extending towards the protostar (e.g. Raga 
et al. 2004). However, in these models the opening angle 
of the cavities near the protostar tends to be extremely 
narrow, in contrast to the wide angle observed in Cep 
E. Furthermore, the emission in the the bow shock is 
distinctly different, predominantly H2 emission, as seen 
in the IRAC images along the jet axis; but in the wide 
angle outflow region the H2 emission is totally absent as 
seen by the high contrast between the emissions in the 4.5 
/im and 8 /im bands in Fig. 4. The lack of H2 emission 
and the prevalence of scattered photospheric light rule 
out the possibility that the wide angle outflow is part 
of the jet driven bow shock. Thus we can constrain the 
extent of any jet precession to be strictly confined to 
the emission region in the 8 jam image as shown by the 
schematic in Fig. 4. 

Whatever is the cause for the wide angle cavity in 
Cep E (by precessing jets or by wide angle low veloc- 
ity outflow), in Spitzer images it is observable only in 
the scattered light. We can clearly distinguish between 
the jet and the wide angle cavity by using the IRAC col- 
ors for scattered light in the wide angle cavity and H2 
emission in the narrow high velocity jets. In Fig. 4a 
the 8 /im emission (dominated by H2 emission in the 
jets and bow shocks) correlates well with the inner nar- 
row jet features in the 4.5 /im image, but not with the 
biconical emission away from the jet. Thus it is clear 
that biconical low surface emission is distinctly different 
from the jet emission and it is coming from the scattered 
light in a wide angle cavity. In Fig. 4a the geometry 
of the wide angle outflow and the narrow high velocity 
jets in the northern lobe are shown schematically by con- 
ical boundaries. The emission in the region outside the 
inner cone clearly represents the scattered light in the 
wide angle outflow and that from inside the inner cone 
boundary represents the jets and their associated bow 
shocks. To study the spectral characteristics of the scat- 
tered photospheric emission from the protostar in the 
wide angle outflow cavity we were able to separate the 
scattered emission from the jet emission by computing 
the flux densities at each IRAC band integrated over the 
area between the inner and outer boundaries of the out- 
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flow cone as shown in Fig. 4a. The scattered hght flux 
density spectral plot is shown in Fig. 4b. This spec- 
trum is distinctly different from those in Fig. 6a for the 
"knots" in the jets and bowshocks. On the other hand, 
between 4.5 /im and 8 /im it is suggestive of the photo- 
spheric (at T=4300K) emission spectrum (shown in Fig. 
9) obtained from the SED analysis (see section 3.4), if we 
take into consideration the effects of higher extinction at 
3.6 /im. The K-band extinctions are in the range of 2.1 
to 2.4 mag across the northern lobe and 1.4 mag for the 
southern lobe (Smith et al. 2003). The lower flux den- 
sity at 3.6 /im, even after the extinction correction (using 
the near-IR to mid-IR extinction relations given by In- 
denbetouw et al. 2005; Chapman et al. 2009), would 
still be consistent with the examples of the SEDs for the 
scattered light which show a dip near 3 /im due to the 
scattering albedo (Whitney et al. 2003b). A compari- 
son between the stellar and the scattered light spectra 
will require detailed modeling including the grain prop- 
erties and the outflow geometry (Whitney et al. 2003a 
& b) as presented for the bright scattered light emission 
from NIR to MIR in the bipolar outflow cavities of the 
Class object L1527 (Tobin et al. 2008). However, it is 
possible that, in addition to the complex inner structure 
generated by the outflow cavities, some of the observed 
scattered light could result from the outer layers of the 
extended disk which forms from infalling matter (Tobin, 
Hartman, & Loinard, 2010). 

3.2.2. H2 Jet Component 

Most of the thermal energy associated with the shocks 
in the outflow is radiated away through the H2 emission 
in the rotational transitions; in the LI 157 outflow, Nisini 
et al. (2010) estimate about 50% of the total shock radi- 
ated energy is in H2 emissions. The similarity between 
the jet/bowshock emissions along the jet axis in all the 
IRAC images is due to the fact that the mid IR jet emis- 
sion is dominated by the H2 pure rotational transitions 
and each IRAC band contains several of the H2 lines. 
(See Neufeld & Yuan (2008) and Ybarra & Lada (2009) 
for the details on their contributions to the band inten- 
sities.) Below we list the H2 lines in each IRAC band 
and their upper state energy in degrees K, as a measure 
of the excitation required for their emission: 
3.6 /im : H2 v = 1 - 0(5)[6952K]; 0(6)[7584K]; 
0(7)[8366K]; H2 v = - S(13)[17445K] 
4.5 /im : H2 v = - S(12)[15542K]; S(11)[13704K]; 
S(10)[11941K]; S(9) [10362K] 

5.8 /im : H2 v = - S(8) [8768K]; S(7)[7197K]; S(6) 
[5830K] 

8.0 /im : H2 V = - S(5) [4587K]; S(4)[3475K] 
To demonstrate how these H2 lines could be used to 
study the outflows. Smith & Rosen (2005) show examples 
of synthetic images of H2 emission in each IRAC band 
using model images obtained in hydrodynamical simula- 
tions of dense supersonic molecular jets. In addition to 
the lines listed above. Smith & Rosen (2005) include in 
each IRAC band a few other H2 v = 1 - 1 and H2 v 
= 2-1 transitions. The 3.6 /im band is dominated by 
emission from vibrationally excited H2 lines with higher 
mean excitation energies, while the other bands are dom- 
inated by rotational lines. 

The H2 emission in a few selected outflows have been 
well studied using the H2 pure rotational lines from S(0) 



to S(7) with IRS mapping (e.g. Neufeld et al. 2006; 
Nisini et al. 2010). No IRS mapping data is available 
for Cep E, however, Noriega- Crespo et al. (2004b) used 
the H2 rotational transition lines in the IRS spectrum 
at a single pointing near the brightest emission in the 
northern lobe (see Fig. 6) to estimate the mean excita- 
tion temperature (686 K) and a number density (6x10^ 
cm~^). Moro-Martm et al. (2001) used the pure rota- 
tional H2 0-0 lines S(2) to S(7) intensities in 7 discrete 
positions observed by ISOCAM CVF and their excita- 
tion diagrams, to model the temperatures and densities 
in the bowshocks. The mid-IR H2 emission is inter- 
preted in terms of shocks produced by the high velocity 
atomic/ionic jets. High-density C-shock models can best 
account for the brightness of the H2 emission (Smith 
et al. 2003; Dionatos et al. 2010). In Cep E by com- 
parison with C-type shock models, the H2 line ratios 
observed in the IRS spectrum give a mean shock veloc- 
ity of 20 km s~^ (Noriega-Crespo et al. 2004b). Since 
the intensities of the H2 lines within each IRAC band 
are sensitive to a wide range of excitation conditions, the 
relative intensities in each band can be used as a diagnos- 
tic of the physical conditions in the shocks. As demon- 
strated by Ybarra & Lada (2009) we can therefore, use 
the IRAC colors: [3.6] - [4.5] versus [4.5] - [5.8], to study 
the shocked H2 at sufficiently high temperature (2000K - 
4000K) and neutral atomic hydrogen density (10^ - 10^ 
cm~^) in jets/outflows. 

The 8 /im image clearly delineates the outer boundary 
of the jets and their associated bowshock emissions as 
shown by the schematic in Fig. 4. In all IRAC bands 
within this region a number of jet features with resolved 
bowshock structures are identified along both north and 
south lobes. Near-IR high resolution H2 1-0 S(l) imag- 
ing has shown that the lobes consist of numerous emis- 
sion knots (Eisloffel et al. 1996; Ladd & Hodapp, 1997). 
Smith et al. (2003) estimate about 20 bow shocks, close 
to paraboloidal in shape, are required by the model to ex- 
plain the integrated emission for each lobe. As indicated 
in the left panels in Figs. 5a & 5b, in the IRAC 8 /im 
image we identify at least 22 such resolved jet/bowshock 
components. As seen in Fig. 1 the emissions in both 
lobes along the jet axis in all four IRAC bands appear to 
be morphologically very similar. Although all the knots 
identified at 8 /im (Fig. 5) have counterparts in all other 
bands, their relative intensities vary widely. Further- 
more, in the 4.5 /im image we could easily identify at 
least an additional seven knots (indicated as Al to A7 in 
Fig. 5) which appear less prominently at 8 /im. To char- 
acterize fully the shock conditions in these II2 knots and 
bowshock features and to bring out their differences, we 
constructed the IRAC channel difference images (Figure 
5): 4.5 /im - 3.6 /im, 5.8 microns - 4.5 /im, and 8 /im - 4.5 
/im. The difference image 4.5 /im - 3.6 /im is more like 
a continuum subtracted H2 line emission (assuming the 
3.6 /im band behaves like a continuum filter with little H2 
line emission) and resembles the 2.12 /im image (Eisloffel 
et al. 1996; Ladd & Hodapp, 1997), tracing the brightest 
emission closer to the leading bowshock edges on small 
scales, while the difference images, 5.8 /im - 4.5 /im and 
8 /im - 4.5 /im, show strong color differences determined 
by the shock excitation along the wings, spatially away 
from their leading edges. 
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To obtain a more quantitative description of the knots 
identified in Fig. 5 we measured their flux densities in 
each IRAC band by integrating over a 1.5^' x 1.5'' area 
and these are also shown as spectral plots in Fig. 6a. 
Though the images have sub-arcsec resolution, to obtain 
a more robust estimation of the flux densities we used 
a larger area of integration. These spectra indicate that 
the excitation conditions vary widely among the knots: 
fourteen (including all seven selected by the 4.5 /im peak 
emission) have apparent spectral peaks at 4.5 /im; ten 
have spectral peaks between 5.8 and 8 /im, while five 
have no apparent peak with fluxes increasing beyond 8 
jam. These differences are also clearly brought out in the 
difference images shown in Fig. 5. The extinction cor- 
rections will increase the fluxes at 3.6 and 4.5 microns 
relative to the others; however, the spectral shapes do 
not change significantly. Our results of such strong 4.5 
/im emission is consistent with the hydrodynamic sim- 
ulations of Smith & Rosen (2005). The knots with the 
highest excitation as inferred from the narrow spectral 
peaks at 4.5 /im (nos. K7, K6 and K2) are all located 
near the brightest tips of the southern and northern lobes 
respectively. While the knots K7 and K6 are among the 
brightest, the knot K2 is only marginally bright. Except 
for these three, there is no other striking correlation be- 
tween the excitation conditions (spectral peak/shape) in 
the knots and their radial distance from the protostar or 
any particular age sequence. 

The IRAC H2 emission spectral data presented here 
will be useful to constrain the C-shock models for each 
jet feature along with those available in the near IR. We 
use the colors [3.6] - [4.5] and [4.5] - [5.8] derived from 
the spectral data for each knot as diagnostics to constrain 
the temperature and densities as discussed by Ybarra & 
Lada (2009). The IRAC color-color plot for the H2 knots 
in Cep E is shown in Fig. 6b. The labels in Fig. 6b rep- 
resent the color-color regime for each H2 knot, identified 
in Fig. 5. The IRAC color-color plot for Cep E knots is 
remarkably similar to those obtained for HH54 (Ybarra 
& Lada, 2009). The majority of the knots, 17 out of 
29, appear to have excitation conditions corresponding 
to temperatures 2500K to 4000K and densities 10^ to 
10^ cm~^. 

The observed bowshock structures (leading edges and 
wings) projected on the sky and the subsequent interpre- 
tation of their shock layers emitting in different IRAC 
bands is critically dependent on the observer's viewing 
geometry with respect to the shock. Taking advantage 
of the sub-arcsec resolution in the deconvolved IRAC im- 
ages we can study in some greater detail the spatial char- 
acteristics in each image to resolve shock geometry. In 
Fig. 7 we show a blowup of the IRAC images near the tip 
of the NE and SW lobes. We use the geometrical shape 
of the 4.5 /im and 8 /im emissions of the shock features 
combined with their relative displacements to "infer" the 
direction the shocks which is likely to be consistent with 
the schematic bowshock structure, for example as shown 
in Fig. 13 in Moro-Martfn et al. (2001). In such cases 
we expect the atomic/ionic emission lines to originate at 
the tip of the bowshock (in the "Mach disk" and stagna- 
tion tip) where the J- shocks are dominant. But, the H2 
line emissions originate farther down along the bowshock 
wings where the shocks are oblique and more likely to be 
C-type. For a better understanding of the relationship 



between the structural complexity in the IRAC emissions 
and their colors with the atomic/ionic jet head, we plot 
in Figs. 7a & 7b the atomic jet traced by the MIPS 24 
/im emission from the [Fe II] and [S I] lines (discussed 
below in section 3.3). The 24 /im emission contour de- 
fines the location of the jet head (Mach disk/ stagnation 
tip) which is driving the bowshock. The high excitation 
H2 lines in the 4.5 /im band are likely to be relatively 
closer to the "tip" than the low excitation traced by 8 
/im band. Schematically, in Figs. 7a & 7b the 4.5 /im in- 
tensity contours overlaid on the 8 /im grey scale delineate 
the bowshock closer to its leading edge (the tip) while the 
8 /im emission traces the projected spatial structure the 
bowshock wings away from the tip. The bright and dark 
bands seen across several knots in the 5.8/im - 4.5/im 
and 8 /im - 4.5/im difference images (eg. nos. 9, 10,12, 
K4, K7, 21 in the right two panels in Figs. 5a & 5b) trace 
the projected separation between the bowshock tips and 
wings. The geometry in two systems of knots consisting 
of nos. 3/5/6 (Fig. 5a and 7a) and K5/18/19 (Fig. 5b) 
respectively, seems consistent with each one describing a 
bowshock tip with wings on either side, moving nearly in 
the plane of the sky. The two outermost arc-like emission 
features in the 4.5/im emission (Fig. 7a) are also detected 
at 2.12 /im (Ladd & Hodapp, 1997). However, the out- 
ermost arc-like emission at 8 /im (Fig. 7a) is located 
between the two 4.5 /im arcs, behind the leading outer 
4.5/im arc and slightly ahead of the second 4.5/im arc, 
with ~ 1.5'' (1000 AU) separation between them. The 
8 /im arc does not seem consistent with bowshock wings 
but rather with a shocked layer between two successive 
leading shocks described by the 4.5/im and 2.12/im arcs 
and originating in two successive ejections. It is inter- 
esting to note that the 4.5 /im outermost arc has a flux 
density spectrum (Kl in Fig. 6a) resembling that of the 
scattered light in the wide angle outflow cavity (Fig. 4b). 
It is, therefore, possible that this outermost arc indeed 
represents scattered light in the cavity produced by the 
NE jet. However, the continuum subtracted 2.12 /im H2 
line emission (Ladd & Hodapp, 1997) strongly favors its 
identifications as a bowshock H2 emission feature. 

The IRAC data for emissions near the SW tip as shown 
in Figs. 5b and 7b suggest a very complex system of 
oblique (C-type) shocks with possible twists in the jet 
kinematics and geometry. The cartoon in Fig. 7b shows 
a plausible scenario for the C-type shocks, seen in pro- 
jection, as traced by IRAC bands. The location of the 
jet head (Mach disk/stagnation tip) which is driving the 
bowshock is indicated by the 24 jam emission contours. 
The multiple C-type shocks suggested by the 4.5 /im and 
8 /im emissions, surrounding the 24 /im jet feature, seem 
to be tracing different pieces of the bowshock and its 
wings, expanding into or out of the plane (but seen in 
projection on the sky), all driven by a single jet-head. 
In Fig. 7b the long arrow marks the general direction 
of the atomic jet while the short arrows indicate oblique 
shocks moving in and out of the plane. We can inter- 
pret the relative outward and inward displacements of 
the 4.5 /im emissions with respect to the 8/im for knots 
K7 and 21, respectively, as due to the viewing geometry 
of the observer and the C-type shock orientations. As 
another example of the geometry of the bowshock with 
respect to the Mach disk, in Fig. 7a we show the asso- 
ciation between the 24 /im atomic/ionic jet feature (#2 
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in Fig. 8a)) and the bowshock traced by knots 4, 5 and 
6 (identified as arc no. 3 in Fig. 7a). 

3.3. Atomic Jets: MIPS 24 fJ^m Compact Sources 

In Fig. 8 we show an overlay of the MIPS 24 fim HiRes 
image (as contours) on the 4.5 /im image. At 24 /im, in 
addition to the protostar, we detect at least three com- 
pact features coincident with the brightest emission in 
the bow shocks as traced by the H2 emission in the 4.5 
/im and other IRAC bands. These compact emissions at 
MIPS 24 /im coincide with the jet heads. Although some 
enhanced emission from the collisionally heated dust in 
the shocks is possible, it is more likely that the 24 /im 
emission is tracing the hottest atomic/ionic gas in the 
bow shocks, possibly similar to that seen in the HH46/47 
system (Velusamy et al. 2007). In contrast to the mor- 
phology of the H2 emission, which are bow shocks with 
extended long arcs, the 24 /im emission is more compact 
and is confined to the tip of the bow shock coinciding 
with the jet head. We can expect to see such differences 
if the 24 /im emission is produced by higher excitation 
(velocity) shocks arising closer to the apex of the jet as 
in the case of J-shocks, while the H2 emission is pro- 
duced further back in the low velocity C-shocks. We 
can, therefore, interpret the 24 /im emission to represent 
that produced in the J-shocks (ionic/atomic shocks) from 
the bright [Fe II] line emissions at 24.51 and 25.98 /im 
within the MIPS pass band, as observed in the case of 
HH46/47. For one of the sources, #1 in Fig. 8a, we can 
use the IRS spectral data to substantiate the fact that 
the 24 /im emission is indeed atomic. In the previously 
published IRS spectrum (Fig. 6 in Noriega- Crespo et al. 
2004b) of this feature we identify only the [S I] ^Pi P2 
line at 24.249 /im, but no [Fe II] line emissions. In order 
to characterize the atomic emissions within the MIPS 24 
/im pass band, we replotted the spectrum using the IRS 
data obtained from the Spitzer Archives. This spectrum, 
shown in Fig. 8b, is consistent with the published spec- 
trum except for the baseline subtraction. The [S I] line 
first detected by Noriega-Crespo et al. (2004b) in this 
position in Cep E is clearly seen in this spectrum. The 
[Fe II] line at 25.98 /im and the H2 S(0) at 28.22 /im 
are also detected with relatively low intensities. We use 
the intensities in this spectrum ([S I] at 0.82Jy and [Fe 
II] at 0.06 Jy) to estimate qualitatively the contribution 
to the 24 /im pass band from the atomic spectral lines 
within the MIPS pass band as indicated in Fig. 6b. The 
total integrated line emission corresponds to 0.64 Jy in a 
single IRS channel (0.17 /im wide which is about 3% of 
the total pass band) appearing at ^ 60% level of the pass 
band response. We would then expect this line emission 
to appear as a 13 mJy continuum source in the MIPS 
24 /im image. Thus, in the HiRes image (with 2'' spa- 
tial resolution) we expect it to have a brightness of ~ 136 
MJy sr~^. This value is in reasonable agreement with the 
observed brightness of 103 MJy sr~^, considering the un- 
certainties in estimating the equivalent intensity in the 
MIPS continuum band for a given channel intensity in 
the IRS spectrum. We therefore conclude that the 24 
/im emission in this jet feature in Cep E is truly atomic 
in origin, primarily from [S I]. Thus after HiRes repro- 
cessing, the Spitzer MIPS 24 /im image is found to be an 
extremely useful diagnostic of the atomic component of 



the protostellar jets. 

For the other two 24 /im sources (#2 and #3 as marked 
by the short arrows in Fig. 8a) seen along the jet emis- 
sion, there are no matching IRS spectra to interpret the 
origin of their emission as atomic/ionic gas in the jet 
heads. However in the case of #3 strong [S II]AA6717/31 
emission is observed by Ayala et al. (2000) in this posi- 
tion which coincides with HH 377. - Furthermore, [S II] 
and [Fe II] have very similar excitation energies (Nisini 
et al. 2002) and therefore, we can expect a significant 
contribution from [S I] and/or [Fe II] to the MIPS 24 
jam emission. Their location along the jet axis coincides 
well with the bright emission knots/jet heads (inferred 
as arising from H2) seen in the IRAC bands. Also their 
positions coincide with the tip of the bow shocks. We 
expect to see [S I] and [Fe II] fine structure line emission 
in the spectra of HH objects, which show the classical 
signatures of collisional excitation similar to those ob- 
served in the IRS spectra of supernova remnants like 
IC 443 (Noriega-Crespo et al. 2009). Therefore, by 
analogy to that identified in the jets of HH46/47 sys- 
tem (Velusamy et al. 2007), we can regard these 24 /im 
emission features to have the same characteristics as the 
source #1 discussed above. Among the three, the source 
#3 in the south lobe is the brightest (peak at 345 MJy 
sr~^) which is also the brightest emission in the IRAC 
bands. The second brightest source #2 (with peak 145 
MJy sr~^) is located in the NE lobe, farther away than 
source #1, but the IRAC emission here is weaker than 
for source #1. In the close vicinity of the sources #1 
and #2 there are multiple bow shocks and jet heads in 
the IRAC bands. However the 24 /im peaks are exclu- 
sively coincident with only one of the jet heads implying 
that not all bow shocks and jet heads have associated 24 
/im sources. In other words, not all jet heads have de- 
tectable atomic/ionic emission. The two outermost bow 
shocks in the northern lobe, which are clearly seen in the 
IRAC bands (see Fig. 7a), do not show any compact 24 
/im sources. Furthermore, these bow shocks do not show 
the bright jet heads as seen in the other cases. 

As pointed out earlier, not all jet heads or knots with 
strong H2 emission in the IRAC bands have associated 
atomic/ionic counterparts. The atomic/ionic emission 
appears to require higher excitation shocks than those 
producing the mid-IR H2 lines. Dionatos et al. (2010) 
discuss in detail the physical conditions required for the 
mid-IR atomic/ionic and H2 line emissions in these jet 
features. They show that the low velocity 10 - 15 km s"-"^ 
shock models that best fit the warm H2 emission are 
unable to reproduce the observed [Fe II] in HH211. It 
would be difficult for us to quantify the differences in H2 
and atomic/ionic emission from the jet heads or knots in 
Cep E without appropriate shock models. However, the 
association between the 24 /im atomic/ionic jet feature 
and the bow shock knots observed in IRAC bands can be 
interpreted in terms of the schematic bowshock model as 
discussed above in section 3.2.2. The 24 /im feature #3 is 
associated with several knots K6, K7, 21 and 22 (in Fig. 
5b and Fig. 7b) and it can be identified as the Mach disk 
driving a system of oblique C-type shocks represented by 
these knots. Similarly, the 24 /im compact source #1 is 
associated with the 4.5 & 8 /im knots no. 10, K3 and 
K4 respectively (Fig. 5a) and can be interpreted as the 
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atomic/ionic jet driving the bow shock represented by 
these knots around it. The 24 /im atomic/ionic jet fea- 
ture #2 brings out even more clearly its association with 
the bow shock and its wings observed in IRAC bands. 
This feature is associated with the knots 4, 5 & 6, and 
the diffuse emission along the bowshock wings (arc no. 
3 as shown in Fig. 7a). 

So far, we have IRS spectral data for only two 24 
/im compact sources associated with the atomic jet head 
(knots), but they seem to have very different atomic spec- 
tral characteristics. The jet head in the HH46/47 system 
shows a strong [Fe II] line while that in Cep E shows a 
strong [S I] line and very weak [Fe II]. Such differences 
suggest that the physical conditions in the jet heads or 
knots can vary widely. The IRS spectral line intensities 
in the jet heads or knots, similar to those in Cep E and 
HH46/47, in HH211 (Dionatos et al. 2010), HH54, HH 
7-11 (Neufeld et al. 2006), BHR71, L1157, L1448, NGC 
2071, and VLA 1623 (Neufeld et al. 2009), show a wide 
variation in the intensity ratio [Fe II] 26 /im to [S I] 25 
/im. Such variations can be interpreted as electron den- 
sity and temperature effects on the ionization of Fe, as 
well as differences in the collisional excitation. Neufeld 
et al. (2009) find evidence for some spatial segregation 
between the [S I] and [Fe II] emitting regions, suggest- 
ing that [Fe II] is tracing the faster, dissociative shocks. 
Thus the differences between the 24 /im features in Cep E 
and HII46/47 are due to their shock conditions, namely 
shocks in Cep E being less dissociative than in IIH46/47. 

The ISO-LWS spectra of both the north and south 
lobes in Cep E detected strong [O I] 63 /im and [C II] 
158/im emissions (Moro-Martm et al. 2001), which are 
the most efficient coolants in low velocity dissociative 
shocks. The low intensity features in the deconvolved 
MIPS 70 /im image (Fig. 2) are very intriguing, espe- 
cially since they lie along the jet /outflow lobes in Cep 
E. Of special interest to us is the fact that the MIPS 
70 /im band includes the [O I] 63 jam line near the end 
of its passband at ^ 75% response. In the ISO LWS 
spectra the [O I] 63 /im line was detected as the bright- 
est line emission with intensities of 9.3 and 12 x 10~^^ 
W cm~^ at the NE and SW lobes respectively. Such 
emission would be observed in the MIPS 70 /im image 
as continuum sources with total flux densities of 0.58 Jy 
and 0.75 Jy respectively for the NE and SW lobes which 
are higher by a factor 2 to 3 than the features in Fig. 
2. These feature are at intensities at a level of 0.1% of 
the peak (protostar) emission. The SW feature is the 
brighter of the two and in the 70 /im Point Response 
Function (PRF) it coincides with an enhanced sidelobe 
in the Airy ring, at ~ 6.3% level. We believe the MIPS 
70 /im data does not warrant such high dynamic range 
(> 1000) in the HiRes processed image. These features 
appear to be residuals of the Airy ring. Furthermore, 
their locations do not seem to correlate with any of the 
features observed in 24 /im or the IRAC bands. 

3.4. Photometric SED Analysis of Cep E Protostellar 
Emission 

We obtained the Spitzer IRAC and MIPS photometry 
using the mosaic images. For the IRAC data we used 
the artifact corrected BCD's (cbcd images) and Imasks 
(bimsk), with the standard overlap and mosaic modules, 
to create the mosaic images for each of the IRAC bands. 



The mosaics were done using a pixratio=2, which results 
in a pixel scale of 0.6^' per mosaic pixel. IRAC photom- 
etry was obtained by using the phot task in IRAF. An 
aperture radius of 4 pixels (2.4^^), sky annulus between 
4 and 8 pixels, was used to obtain the aperture photom- 
etry. For the 3.6 /im, 4.5 /im, 5.8 /im and 8 jam bands 
we used zero point magnitudes of 18.59, 18.09, 17.48 and 
16.7 respectively, and corresponding aperture corrections 
of 1.213, 1.234, 1.379 and 1.584. 

The Spitzer MIPS BCD data in the 24 /im and 70 
/im bands were used to create the basic mosaics using 
the relevant templates in the MOPEX software. Subse- 
quently, photometry was obtained by using the appropri- 
ate default MOPEX/ APEX single- frame templates for 
each band. The recommendations of the Snitzer data 
analysis cookbook recipe 24 were followedjj We note 
that our MIPS photometry agrees with the photometry 
derived by Noriega- Crespo et al. (2004b) to better than 
5%. 

We constructed the spectral energy distribution of the 
Cep E protostar by combining our Spitzer photometry 
with other photometry from the literature. Chini et al. 
(2001) have studied this protostar at millimeter wave- 
lengths, using the fluxes at 450 and 850 /im obtained 
with SCUBA on the JCMT and at 1300 /im with the 
IRAM telescope. The aperture sizes used to derive the 
fluxes are all larger than 10'^ which is larger than the 
aperture used to obtain the MIPS 70 /im photometry. In 
order to model the protostar using the SED fitting tool 
by Robitaille et al. (2007), it is essential to classify the 
available data as "data points" and "upper limits" . The 
data points are used to obtain the best model fitting and 
computing the factor, whereas the upper limits serve 
in constraining the fits. We note that the protostar ap- 
pears as a point source in the Spitzer-IRAC and MIPS 
bands (Fig. 3). Therefore, we use the IRAC and MIPS 
photometry as data points in subsequent modeling. The 
millimeter fluxes from Chini et al. (2001) are obtained 
from apertures larger than 10^' and could be contami- 
nated by the outflow and surrounding cloud. Therefore, 
we use the millimeter data points as upper limits. We 
also use the IRAS 12 /im flux (obtained from an aper- 
ture of 30^' ) as an upper limit. The protostar is not de- 
tected in the 2MASS K band, therefore, an upper limit 
of 16 mag in the K band was used to constrain the SED. 

The SED, characterized by these data points and upper 
limits, is modeled using the SED fitting tool (Robitaille 
et al. 2007). We note that the actual photometry has an 
accuracy closer to l%-2%. However, we assumed pho- 
tometric uncertainties of 10% as a conservative limit to 
minimize the risk of including biases in the modeling. 
The SED fitting tool requires a range of values for the 
distance and interstellar extinction to the source, in or- 
der to scale the models in the grid to the observed data. 
The lower and upper limits to the distance were set at 
0.65 and 0.8 kpc, respectively. The interstellar extinction 
was set to vary between 2 and 8 visual magnitudes. 

The observed SED is well represented by a YSO model 
in the Class O/I evolutionary stage. Fig. 7 shows the 
models fitted to the SED. The detailed physical param- 
eters of the young stellar object system derived from the 
SED fitting are listed in Table 2. The errors quoted in 
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Table 2 correspond to the spread in the parameters of 
the models shown using the grey lines in Fig. 7. Which 
physical parameters are most reliable depends on what 
kind of model best describes the observed SED. This phi- 
losophy is different from other interpretations of YSO's 
based on colors or spectral indices alone (for details see: 
Smith et al. 2010; Gramajo et al, 2010). In a Class O/I 
YSO model, it is the envelope properties, listed in Ta- 
ble 3, that are better constrained, in contrast to the disk 
properties, because the disk is deeply buried inside the 
envelope. The envelope mass for the Cep E protostar is 
~ TMq, with an accretion rate of ^ IO^^Mq yr~^. The 
accreting photosphere embedded inside this envelope is 
currently estimated to have a mass of ~ SMq and a disk 
mass ^ 1O~^M0. Together, these values support the 
conjecture that Cep E will become an intermediate mass 
star of - 4M0 (Froebrich et al. 2003; Froebrich, 2005) at 
the end of star formation. The wide angle flow detected 
in this work may well be the signature of rapid accretion 
and radiation driven wind from the envelope, which is 
more appropriate for intermediate mass stars than low 
mass stars. The inclination angle of the source, obtained 
by model fitting, is ~ 42°zbl7°. Our model fitted value 
agrees, within the uncertainties, with the inclination of 
57° ±7° derived from the radial velocity and proper mo- 
tion (Smith et al. 2003). 

3.5. Infall Envelope as an IR Absorption Band 

Our SED analysis show the existence of an arcmin size 
dense envelope with 7M0 and size 20,000AU (see Table 
3). Such large scale protostellar infall envelopes perpen- 
dicular to the outflow are common in Class protostars 
(cf. Velusamy & Langer 1998; Chang et al. 2010). In 
the Class protostar LI 157, also located in the Cepheus 
flare, an extended ~ 2' flattened envelope perpendicu- 
lar to the outflow was seen in the IRAC band at 8 jam 
in absorption against the background interstellar PAH 
emission (Looney et al. 2007; Tobin et al. 2010). Stutz 
et al. (2008) detected a similar extinction feature (which 
they refer to as a shadow) around the protostar in the 
IRAC 8 /im image of B335 consistent with the flattened 
molecular core. In Fig. 8 we show the mosaic images at 
3.6 /im and 8 /im of a larger region than shown in Fig. 
1. To the south of Cep E we can see a large interstellar 
emission cloud possibly rich in PAH emission, since it is 
not seen in the 3.6 /im image shown for comparison in 
Fig. 8. To the north of this cloud emission, on both 
sides, to the east and the west of the Cep E jet /outflow 
region, we see patches of 8 /im absorption of this PAH 
background. Two patches (indicated by the arrows in 
Fig. 8b) are seen on either side closest to the protostar. 
It seems likely that these are the ends of a continuous 
absorption band which is masked in the middle by the 
protostar and its bright "side lobes" (Airy rings) around 
it. It is possible that these absorption features are sim- 
ilar to those seen in LI 157 by Looney et al. (2007) but 
is seen less clearly here because of the confusion by the 
emission in the sidelobes of the protostar. 

Though the HiRes deconvolved image (Fig. 1) is free 
from the sidelobe contamination from bright sources in 
the image, we cannot use it for tracing the absorption 
feature. The HiRes processed images do not contain 
the background emission to trace this absorption fea- 
tures. Optimal performance HiRes requires that the 



background emission is fully subtracted out in each BCD 
prior to applying the deconvolution. Furthermore, the 
positivity criteria implicit in the deconvolution algorithm 
makes it insensitive to negative intensities. However we 
can still use the HiRes deconvolution on the BCDs with 
the full background emission present and remove the 
sidelobes, but sacrificing the full resolution enhancement. 
For this purpose we stop the HiRes processing with fewer 
than 10 iterations instead of 50 iterations used for opti- 
mal processing. In Fig 8c we show the results of HiRes 
deconvolution after 10 iterations on the data with the 
background present. After the removal of the sidelobes 
the protostellar envelope is traced more clearly by the 
8 /im absorption extending closer to the protostar. We 
interpret this absorption feature as the extended proto- 
stellar envelope. We estimate an angular size of about 
35'', corresponding to 22,000 AU which is in good agree- 
ment with that inferred from the SED analysis. It is 
most likely that this absorption feature is the infall en- 
velope which is feeding the disk and the forming star in 
Cep E. Indeed, Gregersen et al. (2000) have detected the 
infall signature in Cep E using the velocity asymmetries 
observed in HCO+(3-2) with a 26'' beam. A scenario 
with an extended infall region is particularly interesting 
in view of the simultaneous presence of collimated jets 
and wide angle outflow in Cep E. In the recent mod- 
els by Machida et al. (2008) the high velocity jets and 
the low velocity wide angle outflow are driven by two 
different components, namely, the accretion in the com- 
pact protostellar disk driving the jets and the accretion 
from the extended infall envelope driving the wide angle 
outflow. Thus, along with the other examples of B335 
and LI 157, we show that the 8 /im extinction is a useful 
diagnostic of the protostellar envelopes. 

3.6. Simultaneous Jets and Wide Angle Outflows 

In some protostellar objects both wide angle outflows 
and collimated jets are seen, while in others only one is 
observed. These observations raise a number of impor- 
tant questions about jets and outflows, which directly 
impact our understanding of their role in star and disk 
formation. What is the origin of the wide-angle out- 
flows and collimated jets and how do they evolve with 
time? To date, jets and outflows from a protostar have 
various morphological and kinematical properties such 
that they cannot be explained by a single-class model. 
The centrifugal force originates from the rotation of the 
circumstellar accretion disk and the jets are ultimately 
powered by the infall of material, with the rotation and 
magnetic fields playing a crucial role. The magnetocen- 
trifugal origin of jets and their launch from the magne- 
tized accretion disk of the protostar (cf. Ouyed & Pu- 
dritz 1997) are generally accepted, although the detailed 
mechanism is under debate. The wide-angle outflow may 
also be jet driven (cf. Raga & Cabrit 1993, Ostriker et 
al. 2001) or wind driven (cf. Shu et al. 2000). In the 
jet-driven bow shock model, an episodic variation in jet 
velocity produces an internal bow shock driving an in- 
ternal shell, in addition to the terminal shock. Poorly 
collimated flows seen as wide angle outflows can be due 
to extreme precession of the jet (Shepherd et al. 2000) 
and indistinguishable from wind driven wide angle out- 
flows. Other explanations include a turbulent jet model 
(cf. Canto & Raga 1991; Watson et al. 2004), and a 
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circulation model, involving radiation and magnetocen- 
trifugal acceleration and collimation producing heated 
pressure- driven outflows (cf. Lery 2003, Combet et al. 
2006). 

In a unified model for bipolar jets and outflows Shang 
et al. (2006) incorporate the essential features of the pri- 
mary X-wind, which is driven magnetocentrifugally from 
close to the protostar in the interface between its mag- 
netosphere and the associated circumstellar disk. This 
primary wind has an angle- dependent density distribu- 
tion, with a dense axial jet surrounded by a more tenu- 
ous wide-angle wind. The resulting structure shows two 
prominent dense features: a shell of mostly swept-up am- 
bient material and a jet along the axis that is the densest 
part of the primary wind. The shell can be either well- 
coUimated, as observed for the class of jet-like molecular 
outflows, or wide open, as in the classical molecular out- 
flows. The morphology of jets and outflow is shaped 
to a large extent by the ambient mass distribution in 
the collapsing envelope and thus jet-like and wide angle 
outflows are unified only in an evolutionary sequence. 
In the disk-wind model simulations of collapsing, ro- 
tating, magnetized Bonnor-Ebert spheres of molecular 
cores, the highly collimated wind is driven magnetocen- 
trifugally from a wide range of circumstellar disk radii, 
surrounded by a wide-angle wind driven by toroidal mag- 
netic pressure (e.g. Banerjee & Pudritz 2006). In an- 
other unified model Machida et al. (2008) explain the 
simultaneous occurrences of both the jets and wide an- 
gle outflows using the results from their calculation of 
cloud evolution from a molecular cloud core to a stellar 
core, starting with a Bonner-Ebert isothermal core rotat- 
ing in a uniform magnetic field. They find two distinct 
flows, wide angle low velocity and narrow high- velocity 
flows driven by two different components: the collaps- 
ing core and the protostar, respectively. These two dis- 
tinct flows have different degrees of collimation and ve- 
locities: the low-velocity flow (i.e. molecular outflow) 
has a wide opening angle, while the high- velocity flow 
(i.e., optical jet) has a well-coUimated structure. This 
collimation is caused by both the configuration of the 
magnetic field lines around the drivers and their driv- 
ing mechanisms. The low velocity wide-angle outflow, 
mainly driven by the magnetocentrifugal wind mecha- 
nism (disk wind), and guided by hourglass-like field lines; 
and the fast highly- collimated outflow driven by mag- 
netic pressure and guided by straight field lines near the 
protostar. 

The results of the unified model by Machida et al. 
(2008), as seen in their schematic (Fig. 15 in their paper), 
are remarkably similar to the simultaneously present nar- 
row jet and wide angle outflow in Cep E we derive from 
Spitzer observations. This model also appears to be 
viable for the jet/outflow seen in several other objects 
such as in B5-IRS1 in which the wide-angle CO outflow 
(Velusamy and Langer 1998), and the pc scale HH flows 
(Yu et al. 1999) are likely be driven by the extended 
infall and the compact disk, respectively. HH46/HH47 
(Velusamy et al. 2007), Cep A (Torrelles et al. 2011), and 
HH 211 (Hirano et al. 2006) are a few other examples of 
simultaneous jets and wide angle outflows in protostars. 
In outflows of L1448-mm and IRAS 04166+2706, Tafalla 
et al. (2010) observe evidence for the simultaneous pres- 
ence of multiple components consisting of a wide angle 



slow wind (identified by the line-wings) and collimated 
flows (identified by the extremely high velocity (EHV) 
component). A high angular resolution CO(J = 2-1) map 
of the EHV component in 104166 is well explained as re- 
sulting from a series of internal working surfaces travel- 
ing along a collimated jet (S ant iago-G area et al. 2009). 
The large scale morphologies observed in these objects 
are broadly consistent with any of the above models for 
the multiple jet /outflow components. Observations with 
much higher spatial resolution with sufficient sensitivi- 
ties will be needed to probe the launching region of the 
primary winds on a scale of a few AU of the protostar to 
discriminate between various models of jets and outflows. 
Furthermore a large sample of fully mapped jets and out- 
flows will lead to a better understanding of the simulta- 
neously present jets and wide angle outflow sources in 
the context of the evolutionary sequences of protostars. 

4. CONCLUSIONS 

By combining the high sensitivity in the Spitzer im- 
ages and the reprocessing with HiRes we find new jet 
and outflow features in Cep E in addition to those previ- 
ously known; the new results are: (i) wide angle outflow 
seen in the scattered light; (ii) morphological details of at 
least 29 jet driven bow shocks and jet heads or knots; (iii) 
three compact features in 24 /im continuum image iden- 
tified as atomic/ionic line emission coincident with the 
jet heads; (iv) a flattened arcmin size protostellar enve- 
lope seen against the interstellar background PAH emis- 
sion as an absorption band seen across the protostar at 8 
/im. We demonstrate that by separating the protostellar 
photospheric scattered emission in the wide angle cavity 
from the jet emission, we can study directly the scattered 
light spectrum. The simultaneous presence of collimated 
jets and wide angle outflow in Cep E, as shown here, 
is consistent with unified jet-outflow models such as by 
Machida et al. (2008) in which the accretion in the com- 
pact protostellar disk drives the high velocity jets and 
the accretion from the extended infall envelope drives 
the wide angle outflows. We have obtained the H2 emis- 
sion line spectra as observed in all IRAC bands for all 29 
knots and use their IRAC colors as probes of the tem- 
perature and density in the jets and bowshocks. These 
spectra are useful as diagnostics of the C-type shock ex- 
citation of pure rotational transitions of H2 and a few 
H2 vibrational emissions. A shock excitation model tem- 
plate of IRAC intensity spectra taking into account all 
H2 transitions and the spectral response of each IRAC 
band will be a useful resource for interpreting the vast 
Spitzer data base available to study protostellar jets and 
outflows. Detailed modeling of the individual shocks will 
help retrace the history of episodic jet activity and the 
associated accretion on to the protostar. 

Our results for Cep E, which utilize the HiRes algo- 
rithm, show the potential of Spitzer data to study the 
entrained molecular H2 in the high velocity jets as well as 
the wide angle outflow cavities associated with the slower 
winds often traced by CO. Our approach for Cep E builds 
on the earlier work of Velusamy et al. (2007), who iden- 
tified simultaneous jets and wide angle outflows in the 
young protostar HH46/HH47. In addition to the Spitzer 
data, the recent mm and sub-mm spectral line observa- 
tions of outflows and jets from ground (e.g. Tafalla et 
al. 2010; Santiago- Garca et al. 2009) and from Herschel 
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(e.g. Codella et al. 2010; Lefloch et al. 2010), will offer 
new clues to the nature of the EHV gas and its relation 
to the low velocity outflow shells. Though a few selected 
outflows have been well studied in H2 emission using all 
H2 pure rotational lines from S(0) to S(7) with IRS (e.g. 
Neufeld et al. 2006; Nisini et al. 2010) the IRAC im- 
ages available in the Spitzer archives have the potential 
to study the H2 emission in a large sample of outflow 
sources. Quantitative analysis of IRAC colors using sta- 
tistical equilibrium estimates such as those of Ybarra and 
Lada (2009), demonstrate the value of Spitzer images to 
probe the thermal structure of the shocked gas without 
the need of using spectroscopic data. Thus, the Spitzer 
capability to study these features in II2 and scattered 
light, as described here, will be complementary to the 
mm/sub- mm molecular data, thus leading to a better 
understanding of the origin of the atomic jets and the 
fast and slow molecular outflows, and their subsequent 
manifestation in the physical and chemical characteris- 
tics as observed in the resulting shocks. The Spitzer H2 
molecular outflow data is even more relevant now, in the 
light of recent observations (cf. Tafalla et al. 2010) of the 
EHV component (in addition to the wing components) 
in a suite of complex molecules such as CO, SiO, SO, CS, 
HC3N, HCO+, H2CO, HCN, H2O and CH3OH, that can 
be used to model the shock chemistry and excitation con- 
ditions. 

We thank Dirk Froebrich for useful discussions. We 
thank the referee, Alberto Noriega- Crespo, whose critical 
comments helped us to present a more complete and com- 
prehensive picture of Cep E; in particular, for pointing 
out the presence of a second jet and suggesting using the 
IRAC colors. The research described in this paper was 
carried out at the Jet Propulsion Laboratory, California 
Institute of Technology, under a contract with the Na- 
tional Aeronautics and Space Administration. M. S. N. 
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TABLE 1 

Observed Protostar parameters from the HiRes maps 



Band 


Peak 


size 


Flux density 


jim 


MJy sr~i 


// ^ // 


Jy 


3.6 


355 


0.83 X 0.72 


0.0055 


4.5 


1877 


0.73 X 1.00 


0.0358 


5.8 


6465 


0.71 X 0.91 


0.111 


8.0 


10964 


0.78 X 0.85 


0.196 


24 


47597 


1.97 X 2.06 


5.15 


70 


32714 


5.23 X 6.03 


26.4 



TABLE 2 

Results from SED modeling: Protostar 



Age M* R* T* 

log(yr) log(M0) logCRp) logCTp) 
4.33±0.1 0.51±0.01 1.24di0.01 3.64di0.00 



TABLE 3 

Results from SED modeling: Disk and Envelope 





"D env 
^Max 






^disk 


incl 




log(M0yr-i) 


log(AU) 


log(M0) 


log(AU) 


log(M0yr-i) 


o 


log(M0) 


-3.85di0.14 


4.07±0.09 


-2.17di0.73 


1.82ib0.39 


-7.53di0.52 


42±17 


0.87±0.28 
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23^03''18^ 15^ 12^ 09^ 
RA(2000) 



Fig. 1. — Mosaic (lower) and HiRes deconvolved (upper) images in the IRAC bands. In all images the color stretch and the square root 
scaling used are designed to bring out the low surface emissions and thus the brightest emission is saturated. It is readily evident that the 
image enhancement in the HiRes images improves the visualization of the morphologies of jets and outflows by virtue of the increase in the 
resolution to sub-arcsec scale, and by removal of the confusion from the "side-lobes" . 
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Fig. 2. — Mosaic (lower) and HiRes deconvolved (upper) images in the MIPS bands. The arrows in the 24 /im HiRes image indicate 
the features which are distinctly visible only after the HiRes processing and were formerly completely masked by the "side-lobes" in the 
mosaic image. 
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RA(2000) 



Fig. 3. — HiRes maps of the inner 10^^ region around the Cep E protostar. The contour levels are at 0.01%, 0.1%, 1%, 10% and 50% of 
the peak intensity (listed in Table 1) at each band. At 70 fim only the contour at 50% peak intensity is plotted. Note the extremely low 
surface brightness of the outflow-jet features in all IRAC bands, delineated by the low- level contours at intensities below a few percent of 
the peak at the protostar. The crosses in the 3.6 /im map represent the positions of the mm double source (from Moro-Martm et al. 2001). 
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23^03^18^ 16^ 14^ 12^ 10^ 08^ 06^ 

RA(2000) 




A(/^m) 



Fig. 4. — (a) The 8 /im emission (intensity contours) tracing the jet overlaid on the 4.5 /im image which traces both the jet and the wide 
angle outflow in scattered light. The contour intensities are at 1, 2, 4, 8, 16, 32, 64, 128 and 256 MJy sr"-*^ (peak intensity is 10,964 MJy 
sr"-*^). We show a schematic of the bi-conical wide angle outflow (solid lines) and the narrow jets and associated bowshocks (broken lines) 
in the northern lobe. The plausibility of a second H2 molecular jet/outflow (see text) is indicated, (b) Spectral plot of relative intensities 
of the scattered light integrated over the area between the boundaries of the wide angle outflow and the narrow jet as shown in the upper 
panel. 
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a) NE lobe 
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b) SW lobe 
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Fig. 5. — The 8 yum image (left panel), and difference images between IRAC bands (right three panels) showing the spatial variation of 
excitation across the bowshocks inferred from the colors: (a) NE lobe; (b) SW lobe. In H2 a total of 29 jet features, "knots", are identified 
as marked on the images in all panels. The numbered labels 1 to 22 represent the knots identified by the peak emissions in the 8 /xm image. 
The knots identified by Kl to K7 represent those less prominent at 8 yum but are bright at 4.5 fiui. 
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Fig. 6. — (a) Flux density spectral plots of the knots; the labels refer to the position of the knots as marked on Figs. 5a & 5b (see text), 
(b) IRAC color-color plot for the knots identified in Figs. 5a & 5b. The numbered labels identify the knots. The constant temperature 
(solid lines) and density (broken lines), reproduced from Ybarra & Lada (2009), represent statistical equilibrium model calculations of 
IRAC colors for H2 gas in the shocks. 




Fig. 7. — The 4.5 /xm intensity contours (white) overlaid on 8 /xm emission (color image) showing the relative displacement between them. 
The MIPS 24 /im emission (tracing the atomic gas in the jet head/ J- shocks) is shown by selected solid black contours (near peak and half 
power intensities), (a) Tip of NE lobe. Red arcs show a schematic of three bowshocks as inferred from 4.5 /im emissions, (b) Tip of SW 
lobe showing a complex system of oblique C-type shocks. The arcs (dashed lines) represent leading edges of C-type shocks as inferred from 
the 4.5 /im emissions. The short arrows mark the "putative" projected oblique shock directions as "inferred" from the viewing geometry 
and the relative displacement between 8 /xm and 4.5 )U,m intensities (see text). The long arrow marks the direction of the atomic jet. 
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Fig. 8. — (a) 24 /xm emission (contours) overlaid on the 4.5 iim image. The intensity contours are at 4, 8, 16, 32, 64, 128, 256, 512, and 
1024 MJy sr"-*^. The contours near the center represent the protostar with a peak intensity of 47,597 MJy sr"-*^. The long arrow marks 
the emission feature for which the IRS spectrum is shown. The labels 1 to 3 mark the 24 fim emission peaks along the jet. (b) IRS LL 
spectrum of the 24 /im NE emission peak (labeled #1). The atomic/ionic spectral features are marked. The MIPS 24 /im band response 
is also shown. 
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Fig. 9. — SED fitting for the Cep E protostar. The sohd black hne shows the best fitting model. Grey lines show the models that satisfy 
the criteria — X^best < ^- Black dots are data points and triangles are upper limits. The dashed line shows the stellar photosphere of 
the best fit model. 
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Fig. 10. — Extinction at 8 fim through a dense envelope surrounding the protostar. Mosaic images of a large region around the protostar 
showing the interstellar background at 3.6 fim (panel a) and at 8 /im (panel b); the arrows (in panel b) mark extinction near the protostar. 
(c) The 8 /im HiRes deconvolved image after 10 iterations, processed without background subtraction (see text). The ellipse indicates the 
flattened envelope region with high extinction at 8 /xm. 



